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1. Preface 



In this survey I should like to introduce some concepts of algebraic geometry and try to 
demonstrate the fruitful interaction between algebraic geometry and computer algebra 
and, more generally, between mathematics and computer science. One of the aims of 
this article is to show, by means of examples, the usefulness of computer algebra to 
mathematical research. 

Computer algebra itself is a highly diversified discipline with applications to various 
areas of mathematics; we find many of these in numerous research papers, in proceedings 



or in textbooks (cf. 


Buchberger and Winkler (1998) 




Cohen, Cuypers and Sterk (1999), 


Matzat, Greuel and Hiss (1998), 


1SSAC (1988-1998) 


). Here we concentrate mainly on 



Grobner bases and leave aside many other topics of computer al gebra (cf. Davenport, 



Biret and Tournier (1988) , Von zur Gathen and Gerhard (1999)| , Grabmcier, Kaltofen 



and Weispfennig (2000) ). In particular, we do not mention (multivariate) polynomial 
factorisation, another major and important tool in computational algebraic geometry. 
Grobner bases were introduced originally by Buchberger as a computational tool for 
testing solvability of a system of polynomial equations, to count the number of solutions 
(with multiplicities) if this number is finite and, more algebraically, to compute in the 
quotient ring modulo the given polynomials. Since then, Grobner bases have become the 
major computational tool, not only in algebraic geometry. 

The importance of Grobner bases for mathematical research in algebraic geometry 
is obvious and their use needs, nowadays, hardly any justification. Indeed, chapters on 



Grobner bases and Buchberger's algorithm (Buchberger, 1961:) have been incorporated 
in many new textbooks on algebraic geometry such as the books of Cox, Little and 
O'Shea (1992)|, |Cox, Little and O'Shea (1998)| or the recent books of |Eisenbud (1995)| 



Vasconcelos (1998), not to mention textbooks which are devoted exclusively to Grobner 



bases, as Adams and Loustaunou (1994), Becker and Weispfennig (1993), Froberg (1997) 



Computational methods become increasingly important in pure mathematics and the 
above mentioned books have the effect that Grobner bases and their applications become 
a standard part of university courses on algebraic geometry and commutative algebra. 
One of the reasons is that these methods, together with very efficient computers, allow the 
treatment of non-trivial examples and, moreover, are applicable to non-mathematical, 
industrial, technical or economical problems. Another reason is that there is a belief 
that algorithms can contribute to a deeper understanding of a problem. The human idea 
of "understanding" is clearly part of the historical, cultural and technical status of the 
society and nowadays understanding in mathematics requires more and more algorithmic 
treatment and computational mastering. 

On the other hand, it is also obvious that many of the recent deepest achievements in 
algebraic and arithmetic geometry, such as string theory and mirror symmetry (coming 
from physics) or Wiles' proof of Fermat's last theorem, just to mention a few, were neither 
inspired by nor used computer algebra at all. I just mention this in order to stress that no 
computer algebra system can ever replace, in any significant way, mathematical thinking. 

Generally speaking, algorithmic treatment and computational mastering marks not the 
beginning but the end of a development and already requires an advanced theoretical 
understanding. In many cases an algorithm is, however, much more than just a careful 
analysis of known results, it is really a new level of understanding, and an efficient 
implementation is, in addition, usually a highly nontrivial task. Furthermore, having a 
computer algebra system which has such algorithms implemented and which is easy to 
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use, it then becomes a powerful tool for the working mathematician, like a calculator for 
the engineer. 

In this connection I should like to stress that having Buchberger's algorithm for com- 
puting Grobner bases of an ideal is, although indispensable, not much more than having 
+ ,-,*,/ on a calculator. Nowadays there exist efficient implementations of very involved 
and sophisticated algorithms (most of them use Grobner bases in an essential way) al- 
lowing the computation of such things as 

Hilbert polynomials of graded ideals and modules, 
free resolutions of finitely generated modules, 
Ext, Tor and cohomology groups, 

infinitesimal deformations and obstructions of varieties and singularities, 

versal deformations of varieties and singularities, 

primary decomposition of ideals, 

normalisation of affine rings, 

invariant rings of finite and reductive groups, 

Puiseux expansion of plane curve singularities, 



not to mention the standard operations like ideal and radical membership, ideal inter- 
section, ideal quotient and elimin ation of variables. All the above-mentioned a lgorithms 
are implemented in S INGULAR (Greuel, Pfister and Schoncmann, 1990-1998 ), some of 
them also in CoCoA ( |Capani, Nicsi and Robbiano, 1995 ) and Macaulay, resp. Macaulay2 
(Bayer and Stillmann, resp. Grayson and Stillmann), to mention computer algebra sys- 
tems which are designed for use in algebraic geometry and commutative algebra. Even 
general purpose and commercial systems such as Mathematica, Maple, MuPad etc. of- 
fer Grobner bases and, based on this, libraries treating special problems in algebra and 
geometry. 

It is well-acknowledged that Grobner bases and Buchberger's algorithm are responsible 
for the possibility to compute the above objects in affine resp. projective geometry, that 
is, for non-graded resp. graded ideals and modules over polynomial rings. It is, however, 
much less known that standard bases ( "Grobner bases" for not necessarily well-orderings) 
can compute the above objects over the localisation of polynomial rings. This is basically 
due to Mora's modification of Buchberger's algorithm (Mora (1982)) which has been 
modified and extended to arbitrary (mixed) monomial orderings in SINGULAR since 
1990 and was published in |Grassmann ct al (1994)| and in Greuel and Pfister (1996)| . We 
include a brief description in Section 5. 

I shall explain how non well orderings are intrinsically associated with a ring which 
may be, for example, a local ring or a tensor product of a local and a polynomial ring. 
These "mixed rings" are by no means exotic but are necessary for certain algorithms 
which use tag-variables which have to be eliminated later. The extension of Buchberger's 
algorithm to non-well-orderings has important applications to problems in local algebraic 
geometry and singularity theory, such as the computation of 



local multiplicities, 

Milnor and Tjurina numbers, 

syzygies and Hilbert-Samuel functions for local rings 



but also to more advanced algorithms such as 
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classification of singularities, 
semi-universal deformation of singularities, 
computation of moduli spaces, 
monodromy of the GauB-Manin connection. 

Moreover, I demonstrate, by means of examples, how some of the above algorithms 
were used to support mathematical research in a non-trivial manner. These examples 
belong to the main methods of applying computer algebra successfully: 

producing counter examples or giving support to conjectures, 
providing evidence and prompting proofs for new theorems, 
constructing interesting explicit examples. 

The mathematical problems I present were, to a large extent, responsible for the de- 
velopment of SINGULAR, its functionality and speed. 

Finally, I point out some open problems in mathematics and non-mathematical appli- 
cations which are a challenge to computer algebra and where either the knowledge of an 
algorithm or an efficient implementation is highly desirable. 

Acknowledgement: The author was partially supported by the DFG Schwerpunkt 
"Effiziente Algorithmen fiir diskretc Problcmc". Special thanks to C. Losscn and, in 
particular, to T. Kcilcn for preparing the pictures. Finally, I should like to thank the 
referees for useful comments. 

2. Introduction by pictures 

The basic problem of algebraic geometry is to understand the set of points x — 
(xi, . . . , x n ) G K n satisfying a system of equations 

/i(ari,...,a; n ) = 0, 
fk(xi,...,x n ) = 0, 

where if is a field and fi, . . . , fk arc elements of the polynomial ring K[x] = K [x\ , . . . , x n ] . 

The solution set of f± = 0, . . . , fk — is called the algebraic set, or algebraic variety 
of /i , . . . , fk and is denoted by 

v = v(f u ...,f k ). 

It is easy to see, and important to know, that V depends only on the ideal 

k 

I = (h, . . . , f k ) = {/ G K[x] | / = ]T Oi/i, a t G K[x}} 

i=i 

generated by f ly . . . , f k in K[x], that is V = V(I) = {x G K n \ f(x) = V / G /}. 
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The Clebsch Cubic 

This is the unique cubic surface which has §5, 
the symmetric group of 5 letters, as symmetry 
group. It is named after its discoverer Alfred 
Clebsch and has the affine equation 

81(a; 3 + y 3 + z 3 ) 

-I89(x 2 y + x 2 z + xy 2 + xz 2 + y 2 z + yz 2 ) 

+54xyz + 126(xy + xz + yz) 

-9(x 2 + y 2 + z 2 ) - 9(x + y + z) + 1 = 0. 




A Cubic with a ^-Singularity 

Degenerating the Cayley cubic we receive 
a ^-singularity. The affine equation is 

x(x 2 - y 2 ) + z 2 (l + z) + \xy + \yz = 0. 




The Cayley Cubic 

There is a unique cubic surface which 
has four ordinary double points, usually 
called the Cayley cubic after its discoverer, 
Arthur Cayley. It is a degeneration of the 
Clebsch cubic, has §4 as symmetry group, 
and the projective equation is 

ZQZJZ2 + z z 1 z 3 + Z Z 2 Z 3 + Z 1 Z 2 Z 3 = 0. 




The Barth Sextic 

The equation for this sextic was found 
by Wolf Barth. It has 65 ordinary dou- 
ble points, the maximal possible number 
for a sextic. Its affine equation is (with 

(8c + 4)x 2 y 2 z 2 - c 4 (x 4 y 2 + y 4 z 2 + x 2 z 4 ) 
+c 2 (x 2 y 4 + y 2 z 4 + x 4 z 2 ) 

2c -I- 1 2 . 2 1 2 i\2 n 

— (x +y +z -1) =0. 
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An Ordinary Node 

An ordinary node is the most simple sin- 
gularity. It has the local equation 

x 2 + y 2 - z 2 = 0. 




A 5-nodal plane curve of degree 11 
with equation 
-16x 2 + 1048576y 11 - 720896t/ 9 
+ 180224y 7 - 19712j/ 5 + 880y 3 - lly + §, 
a deformation of Aio : y 11 — x 2 = 0. 




Whitney's Umbrella 

The Whitney umbrella is named after Has- 
sler Whitney who studied it in connection 
with the stratification of analytic spaces. It 
has the local equation 




This space curve is given parametrically 
by x — t 4 , y — t 3 , z = t 2 , or implicitly by 

x - z 2 = y 2 - z 3 = 0. 
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Of course, if for some polynomial / 6 K[x], f d \ v = 0, then f\ v = and hence, V = V(I) 
depends only on the radical of J, 

\fl = {/ e K[x] | f e I, for some d}. 
The biggest ideal determined by V is 

/(^) = {/£]f[ a ;]|/( a ;) = 0Vxey}, 

and we have / C \fl C /(V) and = V(y/l) = V(/) = V. 

The important Hilbert Nullstellensatz states that, for K an algebraically closed field, we 
have for any variety V C K n and any ideal J C K[x], 

V = V(J)^>I(V) =VJ 

(the converse implication being trivial). That is, we can recover the ideal J, up to radical, just 
from its zero set and, therefore, for fields like C (but, unfortunately, not for R) geometry and 
algebra are "almost equal" . But almost equal is not equal and we shall have occasion to see that 
the difference between I and \fl has very visible geometric consequences. 

Many of the problems in algebra, in particular, computer algebra, have a geometric origin. 
Therefore, I choose an introduction by means of some pictures of algebraic varieties, some of 
them being used to illustrate subsequent problems. 

The above pictures were not only chosen to illustrate the beauty of algebraic geometric objects 
but also because these varieties have had some prominent influence on the development of 
algebraic geometry and singularity theory. 

The Clebsch cubic itself has been the object of numerous investigations in global algebraic 
geometry, the Cayley and the Dzi-cubic also, but, moreover, since the 7?4-cubic deforms, via 
the Cayley cubic, to the Clebsch cubic, these first three pictures illustrate deformation theory, 
an important branch of (computational) algebraic geometry. 

The ordinary node, also called Ai-singularity (shown as a surface singularity) is the most 
simple singularity in any dimension. The Barth sextic illustrates a basic but very difficult and 
still (in general) unsolved problem: to determine the maximum possible number of singularities 
on a projective variety of given degree. In Section 7.3 we report on recent progress on this 
question for plane curves. 

Whitney's umbrella was, at the beginning of stratification theory, an important example for 
the two Whitney conditions. We use the umbrella in Section 4.2 to illustrate that the algebraic 
concept of normalisation may even lead to a parametrisation of a singular variety, an ultimate 
goal in many contexts, especially for graphical representations. In general, however, such a 
parametrisation is not possible, even not locally, if the variety has dimension bigger than one. 
For curve singularities, on the other hand, the normalisation is always a parametrisation. Indeed, 
computing the normalisation of the ideal given by the implicit equations for the space curve in 
the last picture, we obtain the given parametrisation. Conversely, the equations are derived 
from the parametrisation by eliminating t, where elimination of variables is perhaps the most 
important basic application of Grobner bases. 

Finally, the 5-nodal plane curve illustrates the global existence problem described in Section 
7.2. Moreover, these kind of deformations with the maximal number of nodes play also a promi- 
nent role in the local theory of singularities. For instance, from this real picture we can read off 
the intersection form and, hence, the monodromy of the singularity A 10 by a beautiful theory 
of A'Campo and Gusein-Zade. We shall present a completely different, algebraic algorithm to 
compute the monodromy in Section 6.3. 

For more than a hundred years, the connection between algebra and geometry has turned 
out to be very fruitful and both merged to one of the leading areas in mathematics: algebraic 
geometry. The relationship between both disciplines can be characterised by saying that algebra 
provides rigour while geometry provides intuition. 

In this connection, I place computer algebra on top of rigour, but I should like to stress its 
limited value if it is used without intuition. 
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3. Some problems in algebraic geometry 

In this section I shall formulate some of the basic questions and problems arising in algebraic 
geometry and provide ingredients for certain algorithms. I shall restrict myself to those algo- 
rithms where I am somehow familiar with their implementations and which have turned out to 
be useful in practical applications. 

Let me first recall the most basic but also most important applications of Grobner bases to 
algebraic constructions (called "Grobner basics" by Sturmfels) . Since these can be found in more 
or less any textbook dealing with Grobner bases, I just mention them: 

Ideal (resp. module) membership problem 

Intersection with subrings (elimination of variables) 

Intersection of ideals (resp. submodules) 

Zariski closure of the image of a map 

Solvability of polynomial equations 

Solving polynomial equations 

Radical membership 

Quotient of ideals 

Saturation of ideals 

Kernel of a module homomorphism 

Kernel of a ring homomorphism 

Algebraic relations between polynomials 

Hilbert polynomial of graded ideals and modules 

The next questions and problems lead to algorithms which are slightly more (some of them 
much more) involved. They are, nevertheless, still very basic and quite natural. I should like to 
illustrate them by means of four simple examples, shown in the pictures of this section, referred 
to as Example 1) - 4): 

Assume we are given an ideal I C K[xi, . . . , x n ] by a set of generators fi,...,fh & K[x\. 
Consider the following questions and problems: 



1 Is V(I) irreducible or may it be decomposed into several algebraic varieties? If so, find its 
irreducible components. Algebraically this means to compute a primary decomposition of 
I or of \fl , the latter means to compute the associated prime ideals of I. 

Example 1) is irreducible, Example 2) has two components (one of dimension 2 and one 
of dimension 1), Example 3) has three (one-dimensional) and Example 4) has nine (zero- 
dimensional) components. 

2 Is I a radical ideal (that is, I — vJ) ? If not, compute its radical \fl . 

In Examples 1) - 3) / is radical while in Example 4) \/T = {y 3 — y, x 3 — x), which is much 
simpler than I. In this example the central point corresponds to V((x,y} 2 ) which is a fat 
point, that is, it is a solution of I of multiplicity (= dinijc K[x, y]/(x, y) 2 ) bigger than 
1 (equal to 3). All other points have multiplicity 1, hence the total number of solutions 
(counted with multiplicity) is 11. This is a typical example of the kind Buchberger (resp. 
Grobner) had in mind at the time of writing his thesis. 

3 A natural question to ask is "how independent are the generators fi,...,fk of I?", that 
is, we ask for all relations 

(n, — , rfe) € K[x] k , such that ^ nfi = 0. 



These relations form a submodule of K[x] , which is called the syzygy module of 
fx, . . . , fk and is denoted by syz (7). It is the kernel of the 7<[a;]-linear map 



K[x] k 



K[x]; (n,...,r k )t — y^nfi. 



More generally, we may ask for generators of the kernel of a K[x]-linear map K[x] r — > 
7f[a;] s or, in other words, for solutions of a system of linear equations over K[x]. 
A direct geometric interpretation of syzygies is not so clear, but t here are instance s where 
properties of syzygies have important geometric consequences cf. Schreyer (1986). 
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In Example 1) we have syz (7) = 0, in Example 2), syz (7) = ((— y,x)} C 7f[:r] 2 , in 
Example 3), syz (7) = {{-z,y,0),(-z,0,x)) C K[x] :i and in Example 4), syz (7) C K[x] 4 
is generated by (x, —y, 0, 0), (0, 0, x, —y), (0, x 2 — 1, —y 2 + 1,0). 



Four examples 




• • • 

Example 3) the space curve Example 4) the set of points 

V(xy, xz, yz) V(y 4 - y 2 ,xy A - xy, x s y - xy, x 4 - x 2 ) 

5 A more geometric question is the following. Let V(I') C V(I) be a subvariety. How can 
we describe V(I) \ V(I') ? Algebraically, this amounts to finding generators for the ideal 
quotient 

7 : 7' = {/ G K[x] | fl' C 7}. 
(The same definition applies if 7,7' are submodules of K[x] .) 

Geometrically, V(I : 7') is the smallest variety containing V(I) \ V(I') which is the 
(Zariski) closure of V(I) \ V(I'). 

In Example 2) we have {xz, yz) : (x,y) — z and in Example 3) {xy, xz, yz) : {x,y) = {z, xy) , 
which gives, in both cases, equations for the complement of the z-axis x — y = 0. In 
Example 4) we get 7 : {x,y} 2 = {y(y 2 — l),x(x 2 — 1), (x 2 — l)(y 2 — 1)) which is the zero 
set of the eight points V(I) with centre removed. 

6 Geometrically important is the projection of a variety V(I) C K" into a linear subspace 
K n ~ r . Given generators fi,...,fk of I , we want to find generators for the (closure of 
the) image of V(I) in K n ~ r — {x\xi = • • • = x r = 0}. The image is defined by the ideal 
I fl K[x r +i, ■ ■ ■ ,x n ] and finding generators for this intersection is known as eliminating 
xi,..., x r from fi,... ,f k . 
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Projecting the varieties of Examples 1) - 3) to the (x, ?/)-plane is, in the first two cases, 
surjective and in the third case it gives the two coordinate axes in the (x, i/)-plane. This 
corresponds to the fact that the intersection with K[x, y] of the first two ideals is 0, while 
the third one is xy. 

Projecting the 9 points of Example 4) to the a;-axis we get, by eliminating y, the poly- 
nomial x 2 (x — l)(x + 1), describing the three image points. From a set theoretical point 
of view this is nice, however it is not satisfactory if we wish to count multiplicities. For 
example, the two border points are the image of three points each, hence should appear 
with multiplicity three. That this is not the case can be explained by the fact that elimina- 
tion computes the annihilator ideal of K[x,y]/I considered as ^[xj-module (and not the 
Fitting ideal). This is related to the well-known fact that elimination is not compatible 
with base change. 

7 Another problem is related to the Riemann singularity removable theorem, which states 
that a function on a complex manifold, which is holomorphic and bounded outside a sub- 
variety of codimension 1, is actually holomorphic everywhere. This is well-known for open 
subsets of C, but in higher dimension there exists a second singularity removable theorem, 
which states that a function, which is holomorphic outside a sub-variety of codimension 2 
(no assumption on boundedness) , is holomorphic everywhere. 

For singular complex varieties this is not true in general, but those for which the two 
removable theorems hold are called normal. Moreover, each reduced variety has a nor- 
malisation and there is a morphism with finite fibres from the normalisation to the variety, 
which is an isomorphism outside the singular locus. 

The problem is, given a variety V(I) C K n , find a normal variety V(J) C K m and a 
polynomial map K m — > K n inducing the normalisation map V(J) — > V(I). 
The problem can be reduced to irreducible varieties (but need not be, as we shall see) and 
then the equivalent algebraic problem is to find the normalisation of K[xi, . . . ,x n ]/I, that 
is the integral closure of K[x]/I in the quotient field of K[x]/I and present this ring as an 
affine ring K[yi, . . . , y m ]/J for some m and J. 

For Examples 1) - 4) it can be shown that the normalisation of the first three varieties is 
smooth, the last two are the disjoint union of the (smooth) components. The corresponding 
rings are K[xi, £2], K[xi,X2] © -^[^3], K[xi] © K[x2] © -^"[2:3]. The fourth example has no 
normalisation as it is not reduced. 

A related problem is to find, for a non-normal variety V , an ideal H such that V(H) is 
the non-normal locus of V. The normalisation algorithm described below solves also this 
problem. 

In the examples, the non-normal locus is equal to the singular locus. 

8 The significance of singularities appears not only in the normalisation problem. The 
study of singularities is also called local algebraic geometry and belongs to the basic 
tasks of algebraic geometry. Nowadays, singularity theory is a whole subject on its own. 
A singularity of a variety is a point which has no neighbourhood in which the Jacobian 
matrix of the generators has constant rank. 

In Example 1) the whole t-axis is singular, in the three other examples only the origin. 
One task is to compute generators for the ideal of the singular locus, which is itself a 
variety. This is just done by computing sub- determinants of the Jacobian matrix, if there 
are no components of different dimensions. In general, however, we need, additionally, to 
compute either the equidimensional part and ideal quotients or a primary decomposition. 
In Examples 1) - 4), the singular locus is given by (x,y), (x,y,z), (x,y,z), (x,y) 2 , 
respectively. 

9 Studying a variety V(I), I — (/1, . . . , fk), locally at a singular point, say the origin of K n , 
means studying the ideal IK[x]( x \ generated by I in the local ring 

K \x\(x) = |^ I f,9 G K[x],g (xi,...,iE„)j . 

In this local ring the polynomials g with g(0) / are units and K[x] is a subring of 
K[x] (x) . 

Now all the problems we considered above can be formulated for ideals in K \x]( x ) and 
modules over K[x]^ instead of K[x]. 
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The geometric problems should be interpreted as properties of the variety in a neighbour- 
hood of the origin, or more generally, the given point. 



It should not be surprising that all the above problems have algorithmic and computational 
solutions, which use, at some place, Grobner basis methods. Moreover, algorithms for most 
of these ha ve been implemented quite efficiently in several com puter algebra systems, such as 



CoCoA, cf. Capani, Niesi and Robbiano (19951 
and SINGULAR, cf. 



Macaulav2. cf. |Grayson and Stillmann (1996) 
(1990-1998J, the latter also being able to 



Grcucl, Ptister and Schonemann 
handle, in addition, local questions systematically. 

The most complicated problem is the primary decomposition, the latest achievement is the 
normalisation, both being implemented in SINGULAR. 

At first glance, it seems that computation in the localisation iffa;]/^ requires computation 
with rational functions. It is an important fact that this is not necessary, but that basically the 
same algorithms which were developed for K[x] can be used for /("[a;]^ . This is achieved by the 
choice of a special ordering on the monomials of K[x] where, loosely speaking, the monomials 
of lower degree are considered to be bigger. 

However, such orderings are no longer well-orderings and the classical Buchberger algorithm 
would not terminate. Mora discovered, cf. Mora (1982), that a different normal form algorithm, 
or, equivalently, a different division with remainders, leads to termination. Thus, Buchberger's 
algorithm with Mora's normal form is able to compute in if[a;]( x ) without denominators. 

Several algorithms for K[x] use elimination of (some auxiliary extra) variables. But variables 
to be eliminated have, necessarily, to be well-ordered. Hence, to be able to apply the full power 
of Grobner basis methods also for the local ring K[x]/ X ), we need mixed orders, where the 
monomial orderi ng restricted to some varia bles is not a well-ordering, while restricted to other 
variables it is. In |Greuel and Pfister (1996)| and Grassmann et al (1994), the authors described a 
modification of Mora's normal form, which terminates for mixed ordering, and more generally, 
for any monomial ordering which is compatible with the natural semigroup structure. 



4. Some global algorithms 

Having mentioned some geometric problems, I shall now illustrate two algorithms related to 
these problems: primary decomposition and normalisation. 



4.1. PRIMARY DECOMPOSITION 



Any ideal I C R in a Noetherian ring can be written as / = h qi with qi primary ideals (that 

i — 1 

is, qi ^ R and gf £ qi implies g £ qi or f p £ qi for some p > 0). 

This generalises the unique factorisation (valid in factorial rings) / = /f 1 ■ . . . • f Pr with /; 
irreducible, from elements to ideals. In K[x] we have both, unique factorisation and primary 
decomposition and any algorithm for primary decomposition needs factorisation (because a 
primary decomposition of a principal ideal I = (/) is equivalent to a factorisation of /). 

In contrast to factorisation, primary decomposition is, in general, not unique, even if we 
consider minimal decompositions, that is, the associated primes Pi = ^/qi are all distinct and 
none of the qi can be omitted in the intersection. However, the minimal (or isolated) primes, 
that is, the minimal elements of Ass (J) = {pi, . . . ,p r } with regard to inclusion, are uniquely 
determined. The minimal primes are the only "geometrically visible" primes in the sense that 

V(I) = (J V( Pj ) 

Pj EminAss (J) 

is the decomposition of V(I) into irreducible components. A non-minimal associated prime 
Pi £ minAss (/) is called embedded, because there exists a pj £ minAss (I), pj C Pi- This means 
geometrically V(pi) C V(pj), that is, the irreducible component of V(I) corresponding to pi is 
embedded in some bigger irreducible component. 

As an example we compute the primary decomposition of the ideal / = (x 2 y :i —x 3 yz, y 2 z—xz 2 ) 
in SINGULAR, the output being slightly changed in order to save space. 
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LIB "primdec . lib" ; //calling library for primary decomposition 

ring R = 0, (x,y ,z) ,dp; 
ideal I = x2y3-x3yz,y2z-xz2; 
primdecGTZ(I) ; 

==> [1]: [1]: [2]: [1]: [3]: [1]: 

_[l]=-y2+xz _[l]=z2 _[l]=z 

[2] : _ [2] =y _ [2] =x2 

_[l]=-y2+xz [2]: [2]: 

_[l]=z _[l]=z 

_ [2] =y _ [2] =x 

The result is a list of three pairs of ideals (for each pair, the first ideal is the primary component, 
the second ideal the corresponding prime component). The second prime component [2] : [2] is 
embedded in the first [1] : [2]. The first primary component [1] : [1] is already prime, the other 
two are not. 

Hence, / = (y 2 — xz) n (y, z 2 ) D (x 2 , z) and we obtain: 

V{I) = {y 2 - xz = 0} U {y = z 2 = 0} U {x 2 = z = 0} 

(embedded component) 






Primary decomposition 



All known algorithms for primary decompositions in K[x] are quite involved and use many 
different sub-algorithms from various parts of computer algebra, in particular Grobner bases, 
resp. characteristic sets, and multivariate polynomial factorisation over some (algebraic or tran- 
scendental) extension of the field K. For an efficient implementation which can treat examples 
of interest in algebraic geometry, a lot of extra small additional algorithms hav e to be used. In 
particular one sho uld use "easy" splitting as soon and as often as possible, see Decker, Greuej 
a nd Pfister (1998)f 

in SINGULAR the algorithms of Gianni, Trager and Zacharias (1988)| (wh ich was the first 



practical and gene ral primary decomposition algorithm) , the recent algorithm of Shimoyama and 
Yokovama, and some of the homolog' ica.l algebra algor thms for primary decomposition of 



Eisenbud, Huneke and Vasconcelos (1992) have been implemented. For detailed and improved 
version s of these algorithms, together with extensive comparisons, see Decker, Greuel and Pfister] 



(1998)| . 

Here are some major ingredients for primary decomposition: 

1 Reduction to zero-dimensional primary decomposition (GTZ); 
maximal independent sets; 
ideal quotient, saturation, intersection. 



Computer Algebra and Algebraic Geometry 13 



2 Zero-dimensional primary decomposition (GTZ); 
lexicographical Grobner basis; 
factorisation of multivariate polynomials; 
generic change of variables; 
primitive element computation. 



Related algorithms: 



1 Computation of the radical; 
square-free part of univariate polynomials; 
find (random) regular sequences (EHV). 

2 Computation of the equidimensional part (EHV); 
Ext-annihilators; 

ideal quotients, saturation and intersection. 



To see how homological algebra comes into play, let us compute the equidimensional part of 
V(I), that is, the union of all maximal dim ensional components of V(I), or. algebraica lly, the 



intersection of all minimal primes. Following Eisenbud, Huneke and Vasconcelos (1992) 



calculate the equidimensional part of a variety via Ext-groups: 

If c = codim K [ x j (I), th en the equidimensional part of I is the ann ihilator ideal of the module 



Ext c K ^(K[x]/ 1 , K[x]) by Eisenbud, Huneke and Vasconcelos (1992) 



For example, the equidimensional part of V — {xz — yz — 0} is given by the ideal (z) 
ann (Ext^if [z, y, z]/(xz, yz),K[x, y, z])) . 



Using SINGULAR, we obtain this via: 



xz=yz=0 



LIB "homolog.lib"; 

ring r = 0, (x,y,z) ,dp; 

ideal I = xz, yz; 

module M = Ext_R(l,I); 

quotient (M,f reemodule (nrows (M) ) ) ; 

==> _[1] =z 



*=y=0 

Note that module M = Ext_R(i,I) computes a presentation matrix of Ext ! (i?/7, R). Hence, 
identifying a matrix with its column space in the free module of rank equal to the number of rows, 
Ext, 1 (R/I, R) = R n /M with R n = f reemodule (nrows (M) ) and, therefore, Ann (Ext 1 (R/I, R)) = 
M : R n = quotient (M,f reemodule (nrows (M))). 

Above, we used the procedure Ext_R(-,-) from homolog.lib. Below we show that the Ext 
groups can easily be computed directly in a system which offers free resolutions, respectively 
syzygies, transposition of matrices and presentations of sub-quotients of a free module (modulo in 
SINGULAR). Indeed, the Ext-annihilator can be computed more directly (and faster) without 
computing the Ext group itself: 
Take a free resolution of R/I: 

0< — R/I< — R < — R ni < . 

Then consider the dual sequence: 

— > Hom(i?, R) Hom(7T 1 , R) ■ ■ ■ . 

This leads to: 

Ext* (R/I, R) = Ker(rf i )/Im(d i " 1 ) and Ann (Ext^R/I, R)) = Im (cf" 1 ) : Ker(tf). 
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The corresponding SINGULAR commands are: 



int i = 1 ; 
resolution 
module Im 
module Ker 
module ext 
ideal ann 



= res(I,i+l) ; 

= transpose (L [i] ) ; 

= syz (transpose (L [i+1] )) ; 

= modulo (Ker , Im) ; 

= quotient(Im,Ker) ; 



/ / the Ext group 
//the Ext-annihilator 



Since the resolution can be computed by iterated syzygy computation, this is a beautiful 
example of geometric use of syzygies. However, the algorithm is not at all obvious, but based 
on the non-trivial theorem of Eisenbud, Huneke and Vasconcelos. 

4.2. Normalisation 



Another important algorithm is the normalisation of K[x] I I where J is a radical ideal. It can be 



used a s a step in the primary decomposition, as proposed in Eisenbud, Huneke and Vasconceloi 
(1992), but is also of independent interest. Several algorithms have been proposed, especia lly 



by Keidenberg (1975), Stolzenberg (1968), Gianni and Tra ger (1997) , Vasconcelos (1991) 
adl 



had escaped the computer algebra community, however, that prauert and Remmert (197fJ 



It 

lad 

given a constructive proof for the ideal of the non-normal locus of a complex space. Within 
this proof they provi de a normality criterion which is essentially an algorithm for computing 
the normalisation, cf. De Jong (T998). Again, to make the algorithm effi cient needed some extra 

The Grauert-Remmert 



Decker, Greuel, De Jong and Pfister (1998) 



work which is described in 

algorithm is implemented in SINGULAR and seems to be the only full implementation of the 
normalisation. 



Criterion ( Grauert and Remmert, 1971 ): Let R = K[x]/I with I a radical ideal. Let J be a 
radical ideal containing a non-zero divisor of R such that V(J) contains the non-normal locus 
of V(I). Then R is normal if and only if R = Homn(J, J). 

For J we may take any ideal so that V(J) contains the singularities of V(I). Since normali- 
sation commutes with localisation, we obtain 

Corollary: Ann(Homn(J, J) / R) is an ideal describing the non-normal locus ofV(I). 

Now Horrii{(J, J) is a ring containing R and if R ^ Homij( J, J) = Ri we can continue with 
Ri instead of R and obtain an increasing sequence of rings R C Ri C i?2 C . . . . 

After finitely many steps the sequence becomes stationary (because the normalisation of 
R = K[x]/I is finite over R) and we reach the normalisation of R by the criterion of Grauert 
and Remmert. 

Ingredients for the normalisation (which is a highly recursive algorithm): 

1 Computation of the ideal J of the singular locus of the ideal I; 

2 computation of a non-zero divisor for J; 

3 ring structure on Hom(J, J); 

4 syzygies, normal forms, ideal quotient. 

SINGULAR commands for computation of the normalisation: 

LIB "normal. lib"; 
ring S = 0, (x,y ,z) ,dp; 
ideal I = y2-x2z; 
list nor = normal(I); 
def R = nor [1] ; 
setring R; 
normap ; 

==> normap [1]=T(1) 
==> normap [2] =T(1)*T(2) 
==> normap [3] =T (2) ~2 



(M) h-> (s,st,t 2 ) 
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In the preceding picture, R, the normalisation of S, is just the polynomial ring in two variables 
T(l) and T(2). (The "handle" of Whitney's umbrella is invisible in the parametric picture since 
it requires an imaginary parameter t.) 

In several cases the normalisation of a variety is smooth (for example, the normalisation of 
the discriminant of a versal deformation of an isolated hypersurface singularity) sometimes even 
an affine space. In this case, the normalisation map provides a parametrisation of the variety. 
This is the case for the Whitney umbrella: V = {y 2 — zx 2 — 0}. 

5. Singularities and standard bases 

A (complex) singularity is, by definition, nothing but a complex analytic germ (V, 0) 
together with its analytic local ring R — C{x}/I, where C{:r} is the convergent power series 
ring in x — Xi, . . . , x n . For an arbitrary field K let R — K[[x]]/I for some ideal / in the formal 
power series ring AT[[a;]]. We call (V, 0) = (Speci?, m) or just R a singularity (m denotes the 
maximal ideal of the local ring R) and write K{x) for the convergent and for the formal power 
series ring if the statements hold for both. 

If 7 C K[x] is an ideal with / C (x) = (xi, . . . ,x„) then the singularity of V(I) at £ K n 
is, using the above notation, K{x)/I ■ K(x). However, we may also consider the local ring 
K[x\( x )/I ■ K[x]( x ) with if [a;] (a) the localisation of K[x] at (x), as the singularity of V(I) at 
0. Geometrically, for K = C, the difference is the following: C{x} / IC{x} describes the variety 
V(I) in an arbitrary small neighbourhood of in the Euclidean topology while C[x]( x )/IC[x]^ 
describes V(I) in an arbitrary small neighbourhood of in the (much coarser) Zariski topology. 

At the moment, we can compute efficiently only in ^[a;]^) as we shall explain below. In many 
cases of interest, we are happy since invariants of V(I) at can be computed in K[x]^ as well 
as in K{x). There are, however, others (such as factorisation), which are completely different in 
both rings. 



Isolated Singularities 




Aoo : x 2 - y 2 = Doc : y 2 - zx 2 = 
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(V, 0) is called non— singular or regular or smooth if K{x)/I is isomorphic (as local ring) 
to a power series ring K(yi, . . . , yd), or if ^[a;]^) // is a regular local ring. 

By the implicit function theorem, or by the Jacobian criterion, this is equivalent to the fact 
that / has a system of generators gi, . . . , g n -d such that the Jacobian matrix of gi, . . . , g n -d has 
rank n — d in some neighbourhood of 0. (V, 0) is called an isolated singularity if there is a 
neighbourhood W of such that W fl (V \ {0}) is regular everywhere. 

In order to compute with singularities, we n eed the notion of standard basis which is a 



genera lisation of the notion of Grobner basis, cf. Greuel and Pfister (1996), |Greuel and Pfistci 



(1998) . 

A monomial ordering is a total order on the set of monomials {s a |a £ N n } satisfying 
x a >x ^ x a+ ~< > x 0+1 for all a,/3,j £ N". 

We call a monomial ordering > global (resp. local, resp. mixed) if Xi > 1 for all i (resp. Xi < 1 
for all i, resp. if there exist i,j so that Xi < 1 and Xj > 1). This notion is justified by the 
associated ring to be defined below. Note that > is global if and only if > is a well-ordering 
(which is usually assumed). 

Any / £ K[x] \ {0} can be written uniquely as / = cx a + /', with c £ K \ {0} and a > a' for 
any non-zero term c'x a of /'. We set lm(/) = x a , the leading monomial of / and lc (/) = c, 
the leading coefficient of /. 

For a subset G C K[x] we define the leading ideal of G as 

L(G) = (lm( 5 )| ff £Gx{0}>K Hl 
the ideal generated by the leading monomials in G \ {0}. 

So far, the general case is not different to the case of a well-ordering. However, the following 
definition provides something new for non-global orderings: 

For a monomial ordering > define the multiplicatively closed set 

5> := {u £ K[x] x {0} | lm (u) = 1} 

and the A-algebra 

R := LocK[x] := S^K[x] = | / £ K[x],u £ S>}, 

the localisation (ring of fractions) of K [x] with respect to S>. We call LocA"[a;] also the ring 
associated to K[x] and >. 

Note that K[x] C LocA[a;] C K[x\^ x ) and LocA[x] = K[x] if and only if > is global and 
Loc K[x] = K[x]( x ) if and only if > is local (which justifies the names). 

Let > be a fixed monomial ordering. In order to have a short notation, I write 

R := LocK[x] = S^Klx] 
to denote the localisation of K[x] with respect to >. 

Let I C -R be an ideal. A finite set G C / is called a standard basis of I if and only if 
L(G) = L(I), that is, for any / £ I \ {0} there exists a g £ G satisfying lm (<?)|lm (/). 

If the ordering is a well-ordering, then a standard basis G is called a Grobner basis. In this 
case R = K[x] and, hence, Gc/C K[x]. 

Standard bases can be computed in the sam e way as Gro bner bases except that we need a 



different normal form. This was first notice d by Mora (1982) for local orderings (called tangent 



cone orderings by Mora) and, in general, by Greuel and Flister (1996), Grassmann et al (1994) 



Let Q denote the set of all finite and ordered subsets G C R. A map 
NF -.RxG^R, (/,G)^NF(/[G) 
is called a normal form on R if, for all / and G, 

(i) NF (f\G) + =► lm (NF (/|G)) £ L(G), 

(ii) / - NF (/|G) £ (G)r, the ideal in R generated by G. 
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NF is called a weak normal form if, instead of (ii), only the following condition (ii') holds: 

(ii') for each f £ R and each G G Q there exists a unit u G R, so that uf — NF (f\G) € {G)r. 

Moreover, we need (in particular for computing syzygies) (weak) normal forms with standard 
representation: if G = {gi, . . . ,gk}, we can write 

k 

/-NF(/|G) = £)aifli, a t eR, 
i=i 

such that lm (/ — NF (/ |G)) > lm (a^gi;) for all i, that is, no cancellation of bigger leading terms 
occurs among the a^. 



Indeed, if / and G consist of polynomials, we can compute, in finitely many steps, weak 
normal forms with standard representation such that u and NF(/|G) are polynomials and, 
hence, compute polynomial standard bases which enjoy most of the properties of Grobner bases. 

Once we have a weak normal form with standard representation, the general standard basis 
algorithm may be formalised as follows: 

Standardbasis(G,NF) [arbitrary monomial ordering] 

Input: G a finite and ordered set of polynomials, NF a weak normal form with standard 
representation. 

Output: S a finite set of polynomials which is a standard basis of {G)r. 

- S = G" 

= {(/,<?) I /,<?££}; 

- while (P 0) 

choose (/, g) G P; 

P = Px {(/,<?)}; 

/ l = NF(spoly(/, 5 ) | S); 

if (h ± 0) 

P = PU{(h,f) | feS}; 

S = SU{h}; 

~ return S; 

Here spoly(/, g) — x 1 ~ c " f — ;c(gj a - 7 ~^ff denotes the s— polynomial of / and g where x a = 
lm (/), x 13 = lm (g), 7 = lcm(a, 0). 

The algorithm terminates by Dickson's lemma or by the noetherian property of the polynomial 
ring (and since NF terminates). It is correct by Buchberger's criterion, which generalises to non- 
well-orderings. 

If we use Buchberger's normal form below, in the case of a well-ordering, Standardbasis ist 
just Buchberger's algorithm: 

NFBuCHBERGER(f,G) [ well-ordering ] 

Input: G a finite ordered set of polynomials, / a polynomial. 

Output: h a normal form of / with respect to G with standard representation. 

-h = f; 

- while (ft 7^ and exist g g G so that lm (g) | lm (h)) 

choose any such g; 
h = spoly (h,g); 

- return h; 

For a n algorithm to compute a , weak normal form in the case of an arbitrary ordering, we 
refer to |Greuel and Pfistcr (1996) . 



To illustrate the difference between local and global orderings, we compute the dimension of 
a variety at a point and the (global) dimension of the variety. 

The dimension of the singularity (V, 0), or the dimension of V at 0, is, by definition, the Krull 
dimension of the analytic local ring Ov.o = K{x)/I, which is the same as the Krull dimension 
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of the algebraic local ring K\x]i x )/I in case I = {fi, . . . , fk) is generated by polynomials, which 
follows easily from the theory of dimensions by Hilbert-Samuel series. 

Using this fact, we can compute dim( V, 0) by computing a standard basis of the ideal (fi, ■ ■ ■ , fk) 
generated in Lociffa;] with respect to any local monomial ordering on K[x]. The dimension is 
equal to the dimension of the corresponding monomial ideal (which is a combinatorial problem). 

For example, the dimension of the affine variety V = V(yx — y,zx — z) is 2 but the dimension 
of the singularity (V, 0) (that is, the dimension of V at the point 0) is 1: 









1 








V:y(x-1) 


= z(x - 1) 


= 0, 


dim( V, 0) = 


1, dimV 


= 2 



Using SINGULAR we compute first the global dimension with the degree reverse lexicograph- 
ical ordering denoted by dp and then the local dimension at using the negative degree reverse 
lexicographical ordering denoted by ds. Note that in the local ring K[x, y]{ x , y ) (represented by 
the ordering ds) x — 1 is a unit. 



ring R = 0, (x,y,z) ,dp; 
ideal i = yx-y,zx-z; 
ideal si = groebner(i); 
si; 

==> si[l]=xz-z, 
==> si [2] =xy-y 
dim (si) ; 
==> 2 

ring r = 0, (x,y,z) ,ds; 
ideal i = yx-y,zx-z; 
ideal si = groebner(i); 
si; 

==> si[l]=y 
==> si[2]=z 
dim(si) ; 
==> 1 



//global ring 

//leading ideal of i is <xz,xy> 

//global dimension = dim R/<xz,xy> 
//local ring 

//leading ideal of i is <y,z> 
//local dimension = dim r/<y,z> 



6. Some local algorithms 

I describe here three algorithms which use, in an essential way, standard bases for local rings: 
classification of singularities, deformations and the monodromy. 

6.1. Classification of singularities 

In a tremendous work, V. I. Arnold started, in the late sixties, the classification of hypersurfacc 
singularities up to right equivalence. Here / and g £ K{xi, . . . ,x n ) are called right equivalent 
if they coincide up to analytic coordinate transformation, that is, if there exists a local K— 
algebra automorphism (p of K (x) such that / = <p(g). His work culminated in impressive lists of 
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normal forms of singularities and, moreover, in a determinator for singularities which allows the 
determination of the normal form for a given power series ([AGV, 11.16]). This work of Arnold 
has found numerous applications in various areas of mathematics, including singularity theory, 
algebraic geometry, differential geometry, differential equations, Lie group t heory and theoretical 
physics. The work o f Arnold was continued by C.T.C. Wall and others, cf. [Wall (1983)| , |Greuej 
and Kroning (1990)| . 

Most prominent is the list of ADE or simple or Kleinian singularities, which have appeared in 
surprisingly different areas of mathematic s, and still toda y, new connections of these singularities 
to other areas are being discovered (cf. Glreuel (1992) for a survey). Here is the list of ADE 
singularities (the names come from their relation to the simple Lie groups of type A, D and E). 



D k 
E 6 
E 7 

Ek 



x\ +1 + x\ + x\ + ■ ■ 
x 1 {x l l- 2 +x%) + x\ + 



+ x 2 



x l + x 2 + x 3 + ' ' ' + x ni 

X2(xl + x\) + xl H + x 2 

x\ + xi + xi H Vx\. 



k > 1 

k > 4 




A3 -Singularity 



D6-Singularity 



i57-Singularity 



Arnold introduced the concept of "modality", related to Riemann's idea of moduli, into sin- 
gularity theory and classified all singularities of modality < 2 (and also of Milnor number < 16). 
The ADE singularities are just the singularities of modality 0. Singularities of modality 1 are 
the three parabolic singularities: 



E 6 = P S = T33 



x 3 +y a + z 3 + axyz, 
x 4 + y 4 + ax 2 y 2 , 
x 3 + y 5 + ax 2 y 2 , 



Es — J10 = T236 
the 3-indexed series of hyperbolic singularities 

T pqr : x p + y q + z T + axyz, a / 



aV4, 

4a 3 + 27/0, 



i + i<l 
q r 



and 14 exceptional families, cf. Arnold, Gusein-Zade and Varchenko (1985). 

The proof of Arnold for his det erminator is. to a great part constructive, and has been partly 



implemented in SINGULAR, cf. Kriiger (1997). Although the whole theory and the proofs 



deal with power series, everything can be reduced to polynomial computation since we deal 
with isolated singularities, which are finitely determined. That is, for an isolated singularity 
/, there exists an integer k such that / and g are right equivalent if their Taylor expansion 
coincides up to order k. Therefore, knowing the determinacy k of /, we can replace / by its 
Taylor polynomial up to order k. 

The determinacy can be estimated as the minimal k such that 

m fe+1 C m 2 jacob(/) 

where m C K(x\, . . . , x n ) is the maximal ideal and jacob(/) = {df/dxi, . . . , df/dx n ). Hence, 
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this k can be computed by computing a standard basis of m 2 jacob(/) and normal forms of 
m l with respect to this standard basis for increasing i, using a local monomial ordering. How- 
ever, there is a much faster way to compute the determinacy dir ectly from a standard basi s of 



m jacob(/), which is basically the "highest corner" described in Greuel and Pfister (1996). 

An important initial step in Arnold's classification is the generalised Morse lemma, or splitting 
lemma, which says that / o tp(xi, . . . , x„) = x\ + ■ ■ ■ + a; 2 . + g(x r +i, . . . , x„) for some analytic 
coordinate change ip and some power series g £ m 3 if the rank of the Hessian matrix of / at 
is r. 

The determinacy allows the computation of <p up to sufficiently high order and a polynomial g 
as in the theorem. This has been implemented in SINGULAR and is a cornerstone in classifying 
hypersurface singularities. 

In the following example we use SINGULAR to get the singularity T^j^u from a database 
A_L ("Arnold's list"), make some coordinate change and determine then the normal form of 
the complicated polynomial after coordinate change. 



LIB "classify. lib" ; 
ring r = 0, (x,y,z) ,ds; 
poly f = A_L("T[5,7,11] ") ; 
f ; 

==> xyz+x5+y7+zll 

map phi = r, x+z,y-y2,z-x; 

poly g = phi(f) ; 

g; 

==> -X2y+yz2+x2y2-y2z2+x5+5x4z+10x3z2+10x2z3+5xz4+z5+y7-7y8+21y9-35yl0 
==> -Xll+35yll+llxl0z-55x9z2+165x8z3-330x7z4+462x6z5-462x5z6+330x4z7 
==> -165x3z8+55x2z9-llxzl0+zll-21yl2+7yl3-yl4 
classify(g) ; 

==> The singularity ... is R-equivalent to T [p,q,r] =T [5 ,7 , 11] 



Ingredients for the classification of singularities: 



1 standard bases for local and global orderings; 

2 computation of invariants (Milnor number, determinacy, . . . ); 

3 generalised Morse lemma; 

4 syzygies for local orderings. 



Beyond classification by normal forms, the construction of moduli spaces for singularities, 
for varieties or for vector bundles is a pretentious goal, theoretically as well as computatio nal. 
First steps towards th is goal for singularities have been undertaken in T. Bayer (2000) and 



Friihbis-Kruger (2000) 



6.2. DEFORMATIONS 



Consider a singularity (V, 0) given by power series fi,. ■ ■ , fk £ K{xi, . . . , x n ). The idea of de- 
formation theory is to perturb the defining functions, that is to consider power series Fi(t, x), . . . , 
Fk(t,x) with Fi(0,x) = fi(x), where t G S may be considered as a small parameter of a param- 
eter space S (containing 0). 

For t 6 S the power series fi,t(x) = Fi(t,x) define a singularity Vt, which is a perturbation 
of V = Vo for t / close to 0. It may be hoped that Vt is simpler than Vo but still contains 
enough information about Vo- For this hope to be fulfilled, it is, however, necessary to restrict 
the possible perturbations of the equations to flat perturbations, which are called deformations. 

Grothendieck's criterion of flatness states that the perturbation given by the Fi is flat if and 
only if any relation between the fi, say 
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lifts to a relation 

^2Ri(t,x)Fi{t,x) = 0, 

with Ri(x, 0) = Ti(x). Equivalently, for any generator (ri, . - . , rt) of syz (/i, . . . , /&) there exists 
an element (Ri, . ■ . , Rk) € syz (i<i, . . . , Fk) satisfying Ri(0, x) — ri{x). Hence, syzygies with 
respect to local orderings come into play. 

There exists the notion of a semi-universal deformation of (V, 0) which contains essentially 
all information about all deformations of (V, 0). 

For an isolated hypersurface singularity f(xi, ■ ■ ■ , x n ) the semi-universal deformation is given 

by 

T 

3=1 

where 1 =: g\, g%, . . . , g T represent a .ft'-basis of the Tjurina algebra 

K(x)/(f,df/dx 1 ,...,df/dx n ), 

t = dim_R- K{x)/{f, df/dxi, . . . , df/dx n ) being the Tjurina number. 

To compute gi, ■ ■ ■ , g T we only need to compute a standard basis of the ideal (/, , . . . , -§^~) 
with respect to a local ordering and then compute a basis of K[x] modulo the leading monomials 
of the standard basis. For complete intersections we have similar formulas. 




Deformation of E-j in AAi 



For non-hypersurface singularities, the semi-universal deformation is much more complicated 
and up to now no finite algorithm is known in general. However, there exists an algo rithm to 
compute this deformation up to arbitrary high order cf. Laudal (1979), Martin (1998), which is 
implemented in SINGULAR. 

As an example we calculate the base space of the semi-universal deformation of the normal 
surface singularity, being the cone over the rational normal curve C of degree 4, parametrised 
by tn» (t,t 2 ,t 3 ,t 4 ). 

Homogeneous equations for the cone over C are given by the 2 x 2-minors of the matrix: 



x y z u \ 

ni = | e Mat 2x4(K\x, y, z, u, v\). 

y z u v ,' 



SINGULAR commands for computing the semi-universal deformation: 
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LIB "deform. lib"; 

ring r = 0, (x,y,z,u,v) ,ds; 

matrix m[2] [4] = x,y,z,u,y,z,u,v; 
ideal f = minor(m,2); 

versal(f ) ; 
setring Px; 
Fs; 

==> Fs[l,l]=-u2+zv+Bu+Dv 
==> Fs [ 1 , 2] =-zu+y v-Au+Du 
==> Fs[l,3]=-yu+xv+Cu+Dz 
==> Fs[l,4]=z2-yu+Az+By 
==> Fs [ 1 , 5] =yz-xu+Bx-Cz 
==> Fs [ 1 , 6] =-y 2+xz+ Ax+Cy 
Js; 

==> Js[l,l]=BD 
==> Js[l,2]=AD-D2 
==> Js[l,3]=-CD 

The ideal Js = {BD, AD - D 2 ,-CD) C K[A, B, C, D] defines the required base space which 
consists of a 3-dimensional component (D — 0) and a transversal 1-dimensional component 
(B = C = A — D = 0). This was the first example, found by Pinkham, of a base space of a 
normal surface having several components of different dimensions. 

The full versal deformation is given by the map (Fs and Js as above) 

K\[A, B, C, D]]/J s — * K[[A, B, C, D, x, y, z, u, v]]/J s + F s . 

Although, in general, the equations for the versal deformation are formal power series, in 
many cases of interest (as in the example above) the algorithm terminates and the resulting 
ideals are polynomial. 

Ingredients for the semi-universal deformation algorithm: 

1 Computation of standard bases, normal forms and resolutions for local orderings; 

2 computation of Ext groups (cf. 4.1) for computing infinitesimal deformations and obstruc- 
tions; 

3 computation of Massey products for determining obstructions to lift, recursively, infinites- 
imal deformations of a given order to higher order; 

4 one of the main difficulties in point 3 is the necessity to compute a completely reduced 
normal form with respect to a local ordering. In general, such a normal form exists only 
as formal power series. In the present situation, however, the reduction has to be carried 
out only for a subset of the variables in a fixed degree and, hence, the complete reduction 
is finite. 



6.3. THE MONODROMY 

Let / G C{xi, . . . ,x„} be a convergent power series (in practice a polynomial) with isolated 
singularity at and /i = dime C{x} / {f xi , . . . , f Xn ) the Milnor number of /. 

Then / defines in an e-ball B e around a holomorphic function to C, / : B s — > C. 

The simple, counterclockwise path 7 in C around induces a C°°-diffeomorphism of Xt (t 7^ 0) 
(as indicated in the figure) and an automorphism of the singular cohomology group H n (Xt, C) 
which is, by a theorem of Milnor, a /x-dimensional C-vector space. This automorphism 

T : H n (X t ,C) -=* H n (X t ,C) 

is called the local Picard— Lefschetz monodromy of /. We address the problem of computing 
the Jordan normal form of T. 



B=C=0 
A-D=0 
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x = s e n/" 1 (o) 




x t = B E nr 1 (t) 




. — _ 7 ^\ = path around 

•o y* 



The first important theorem is the 
Monodromy theorem (Deligne 1970, Brieskorn 1971): The eigenvalues ofT are roots of unity, 
that is, we have 

where M is a complex matrix with eigenvalues in Q. 

Hence, we are left with the problem of computing the Jordan normal form of M. 

It is not at all clear that the purely topological definition of T allows an algebraic and com- 
putable interpretation. The first hint in this direction is that we can compute dime H n (X t ,C), 
according to Milnor's theorem, algebraically by the formula for [i given above. 

Since X t is a complex Stein manifold, its complex cohomology can be computed, via the 
holomorphic de Rham theorem, with holomorphic differential forms, which is the starting point 
for computing the monodromy. 



To cut a long story short, we just mention, cf. Brieskorn (1970), Greuel (1975) 

H' = n n /df A ff 1 " 1 +dfi 11 - 1 , 

-1 



that 



H" 



Sl n+1 /df Adj- 



ure free C{£}- modules (via /* : C{t} — > C{/} C C{x}) of rank equal to /i. Here (fi*, d) denotes 
the complex of holomorphic differential forms in (C n , 0). H' and H" are called Brieskorn lattices. 
We define the local Gaufi- Manin connection of / as 

V : df AH' = df A Q."/df A d^ 1 > H" , 

X/[df Atu] — [du]. 

Note that \/(df A H') gL df AH', that is, v has a pole at 0. Tensoring with C (t), the quotient 
field of C{t}, we can extend v to a meromorphic connection 



® C(t) 

C{t} 



H' 



® C(i) 

C{t} 



(since df AH' ® C(t) = H" ® C(t)) using the Leibnitz rule v( w J/) : 
With respect to a basis ui, . . . ,o; M of H" we have v( w = J2 a i 

j 



\/{uj)y + ujdy/dt. 
Oj and, for any u> = Ys^iVi, 
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V( tJ ) — ~}2 a jiyi + ~}2^idyi/dt. Hence, the kernel of Vi together with a basis of H" , is the same 

i,i i 

as the solutions of the system of rank \i of ordinary differential equations 

^=-Ay, A = ( a!J ) G Mat(,u x p,C(t)) 
in a neighbourhood of in C. The connection matrix, A = Ait', Ai G Mat (fi x fj,, C), has 

z> — p 

a pole at f = and is holomorphic for t ^ 0. If 0t = (<^i, • • • , <?V) i s a fundamental system of 
solutions at a point t 7^ 0, then the analytic continuation of <j>t along the path 7 transforms <f>t 
into another fundamental system <f>' t which satisfies <j>' t = T v t for some matrix T v G GL(/x, C). 



Fundamental fact (Brieskorn, 1971): The Picard-Lefschetz monodromy T coincides with the 
monodromy T v of the Gaufl-Manin connection. 



Brieskorn used this fact to describe in Brieskorn (1970) the essential steps for an algorithm to 
compute the characteristic polynomial of T. Results of Gerar d and Levelt allowed the extension 
of this algorithm to compute the Jordan normal form of T, cf. Serard and Levelt (1973) . An early 
implementation by Nacken in Maple was not very efficient. Recently, Schulze (1999) implemented 
an improved version in SINGULAR which is able to compute interesting examples. 

The algorithm uses another basic theorem, the 



Regularity Theorem (Brieskorn, 1971): The Gaufi-Manin connection has a regular singular 
point at 0, that is, there exists a basis of some lattice in H" <8> C(t) such that the connection 
matrix A has pole of order 1. 

Basically, if A = A-it^ 1 + Ao + Ait + ■ ■ ■ has a simple pole, then T — e 2n 'A-i is the 
monodromy (this holds if the eigenvalues of A-\ do not differ by integers which can be achieved 
algorithmically) . 

SINGULAR example: 



LIB "mondromy . lib" ; 
ring R = 0,(x,y),ds; 

poly f = x2y2+x6+y6; //example of A'Campo (monodromy is not diagonalisable) 
matrix M = monodromy (f) ; 
print (jordanf orm(M) ) ; 





1/2,1, 


0, 


0, 


0, 


0, 


0,0,0,0, 


0, 


0, 


0, 




0, 


1/2,0, 


0, 


0, 


0, 


0,0,0,0, 


0, 


0, 


0, 




0, 


0, 


2/3,0, 


0, 


0, 


0,0,0,0, 


0, 


0, 


0, 




0, 


0, 


0, 


2/3,0, 


0, 


0,0,0,0, 


0, 


0, 


0, 




0, 


0, 


0, 


0, 


5/6,0, 


0,0,0,0, 


0, 


0, 


0, 




0, 


0, 


0, 


0, 


0, 


5/6,0,0,0,0, 


0, 


0, 


0, 




0, 


0, 


0, 


0, 


0, 


0, 


1,0,0,0, 


0, 


0, 


0, 




0, 


0, 


0, 


0, 


0, 


0, 


0,1,0,0, 


0, 


0, 


0, 


==> 


0, 


0, 


0, 


0, 


0, 


0, 


0,0,1,0, 


0, 


0, 


0, 


==> 


0, 


0, 


0, 


0, 


0, 


0, 


0,0,0,7/6,0, 


0, 


0, 


==> 


0, 


0, 


0, 


0, 


0, 


0, 


0,0,0,0, 


7/6,0, 


0, 


==> 


0, 


0, 


0, 


0, 


0, 


0, 


0,0,0,0, 


0, 


4/3,0, 




0, 


0, 


0, 


0, 


0, 


0, 


0,0,0,0, 


0, 


0, 


4/3 



Ingredients for the monodromy algorithm: 



1 Computation of standard bases and normal forms for local orderings; 

2 computation of Milnor number; 

3 Taylor expansion of units in K[x]^ up to sufficiently high order; 

4 computation of the connection matrix on increasing lattices in H" Cg> C(£) up to sufficiently 
high order (until saturation) by linear algebra over Q ; 
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5 computation of the transformation matrix to a simple pole by linear algebra over Q ; 

6 factorisation of univariate polynomials (for Jordan normal form). 

The most expensive parts are certain normal form computations for a local ordering and the 
linear algebra part because here one has to deal iteratively with matrices with several thousand 
rows and columns. It turned out that the SINGULAR implementation of modules (considered 
as sparse matrices) and the Buchberger inter-reduction is sufficiently efficient (though not the 
best possible) for such tasks. 



7. Computer algebra solutions to singularity problems 

We present three examples which demonstrate, in a somewhat typical way, the use of computer 
algebra as stated in the preface: 

1 producing counter examples; 

2 providing evidence and prompting proofs for new theorems; 

3 constructing interesting explicit examples; 

7.1. EXACTNESS OF THE POINCARE COMPLEX 



The first application is a counterexample to a conjectured generalisation of a theorem of 
Baito (1971) which says that, for an isolated hypersurface singularity, the exactness of the 
Poincare complex implies that the defining polynomial is, after some analytic coordinate change, 
weighted homogeneous. 
Theorem ( JSaito, 197l| ): 

If f : C — ► C has an isolated singularity at 0, then the following are equivalent: 



1 X = f 1 (0) is weighted homogeneous for a suitable choice of coordinates. 

(df_ 

\dx t/ 



2 fi = t where fj, = dime C{x} / 
df 



is the Milnor number and 



t = dime C{x}/ f 



dxi 



the Tjurina number. 



3 The holomorphic Poincare complex 

— > C — > Ox «x 

is exact. 



A natural problem is whether the theorem holds also for complete intersections X = / 1 (0) 
with / = (/i,..., f k ) :C n+k - ■ 

T, 

M = X>1) 1_1 dim c CI*}/)/: 



c . Again we have a Milnor number fi and a Tjurina number 

/. 



Ji- 



d(x jl , 



: dim c C{x}V(/i, . . . , f k )C{x} k + Df(C{x} n+k )) 



Theoretical reduction (Greuel, Martin and Pfister, 1985): 

If X is a complete intersection of dimension 1, then (1) <4> (2) (3). 
If k — 2, then (3) =>■ (2) if fi = dime fix — dime fix and i//i,/a are weighted homogeneous. 

1 PflsXeT arid ScliOliemami (1989)' showed that (3) => (2) does not hold in general: 

fi=xy + z ~ , /2 = xz + y k ~ x + yz' 2 (4 < I < k, k > 5) 
is a counterexample. 

The proof uses an implementation of the standard basis algorithm in a forerunner of SINGU- 
LAR and goes as follows: 
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1 Compute /x, dime &x, dime &x to show that Q.' x is exact; 

2 compute r. 

One obtains y, = t + 1, that is, X is not weighted homogeneous. 

To do this we must be able to compute standard bases of modules over local rings. 

The countere xample was found through a computer search in a list of singularities classified 



by Wall (1983) 



7.2. Zariski's multiplicity conjecture 

The attempt to find a counterexample to Zariski's multiplicity conjecture — which says that 
the multiplicity (lowest degree) of a power series is an invariant of the embedded topological type 
— led, finally after many experiments and computations, to a partial proof of this conjecture. 
For this, an extremely fast standard basis computation for O-dimensional ideals in a local ring 
was necessary. 



The following question was posed by Zariski (1971) in his retiring address to the AMS in 
1971. 

Let f = ^2 c aX a G C{x\, . . . , x n }, /(0) = 0, be a hypersurface singularity, and let mult (/) := 
min{|a| | c a ^ 0} be the multiplicity. 

We say that / and g are topological equivalent, / ~ g, if there is a homeomorphism 

(B,f 1 (D)nB,0) (B lS - 1 (0)nB,0) 




Zariski's conjecture may be stated as: / ^ g mult (/) = mult (g). 

The result i s known to be tr ue for curves (Zariski, Le) and weighted homogeneous singularities 



Greuel (1986), O'Shea (1987) 



Our attempt to find a counterexample was as follows: 
Consider deformations of / = fo: 

ft(x) = f(x) +tg(x,t), \t\ small. 

Then use the theoretical fact proved by Le and Ramanujam: 

fo ~ p ft => K/o) = K/0 

("<=" holds also, except for n — 3, where the answer is still unknown) where /i(/o) respectively 
H(ft) are the Milnor numbers. 

We tried to construct a deformation f t of fo where the multiplicity mult(/t) drops but the 
Milnor number /u(/t) is constant. 

Our candidates (a, b, c 6 N) came from a heuristical investigation of the Newton diagram, one 
being the following series: 

/a , b . 3c c+2 c— 1 . c— 1 e— 1 3 , c— 2 c/ 2 . , \2 i ^ t\t 
t — x + y + z + x y + x y z + x y (y + tx) , a, o, c G N. 



Computer Algebra and Algebraic Geometry 27 



Obviously, the multiplicity drops. Computing fi with SINGULAR, we obtain for (a, b, c) = 
(37,27,6): /u(/o) = 4840, i-i(ft) = 4834, thus fo and ft are (unfortunately) not topologically 
equivalent. 

Since the Milnor numbers of possible counter examples have to be very big, we need an 
extremely efficient implem entation of standard bases. For this, the "highest corner" method of 
Greuel and Pfister (1996) was essential. 

Trying many other classes of examples, we did not succeed in finding a counter example. 
However, an analysis of the examples led to the following 
Partial proof of Zariski's conjecture (Greuel and Pfister, 1996): 
Zariski 's conjecture is true for deformations of the form 

ft = gt(x,y) + z 2 h t (x, y), mult (g t ) < mult (/ ). 

There is also an invariant characterisation of the deformations of the above kind. The general 
conjecture is, up to today, still open. 



7.3. CURVES WITH MAXIMAL NUMBER OF SINGULARITIES 



Let C C be an irreducible projective curve of degree d and f(x, y) = a local equation 
for the germ (C, z). Let n{C, z) — dime C{x, y}/(f x , f y ) be the Milnor number of C at z. 

Since the genus of C, g(C) = ( d-1 H d ~ 2 ) _ 5(C) is non-negative (where 8(C) = 8(C,z), 

zee 

8(C,z) — dime R/R, R ~ C{x,y} / (/) and R the normalisation of R), C can have, at most, 
(d — l)(d — 2)/2 singularities. 

It is a classical and interesting problem, which is still in the centre of theoretical research, 
to study the variety V = Vd{Si, . . . , SV) of (irreducible) curves C C P c of de gree d having 
exactly r singularities of prescribed (topological or analytical) type Si, . . . , SV. Among the most 
important questions are: 

Is V 7^ (existence problem)? 

Is V irreducible (irreducibility problem)? 

Is V smooth of expected dimension (T-smoothness problem)? 

A complete answer is only known for nodal curves, that is, for Vd(r) — Vd{Si, . . . , S r ) with Si 
ordinary nodes (Ai-singularities): 

Severi (1921): V d {r) / and T-smooth <^ r < (d ~ 1) 2 (d ~ 2) . 

Harris (1985): V d (r) is irreducible (if / 0). 

Even for cuspidal curves a sufficient and necessary answer to anv of the above auestions is 
unknown. 



X 



A 4- nodal plane curve of degree 5, with 
equation x 5 ~ fa; 3 + ^x - |y 3 + ±y = 0, 
which is a deformation of Eg : x J — y 3 = 0. 




A plane curve of degree 5 with 5 cusps, 
the maximal possible number. It has the 
equation ±f-x 4 y - § §-x 2 y 3 + gt/ 5 - 20x 4 - 
V + f V 4 ~ 12x 2 y + f y 3 + 32x 2 = 0. 
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Concerning arbitrary (topological types of) singularities, we have the following existence the- 
orem, which is, with respect to the exponent of d, asymptotically optimal. 



Theorem: ( |Greuel, Lossen and Shustin, 1998 ; Lossen, 1999). 
" (d + 2Y 

Vd(Si, . . . , S r ) 7^ l /X][=i ^{Si) ^ — 46~~ an ^ a dditi ona l conditions for the five "worst" 
singularities. 

In case of only one singularity we have the slightly better sufficient condition for existence, 

The theorem is just an existence statement, the proof gives no hint how to produce any 
equation. Having a method for constructing curves of low degree with many singularities, Lossen 
was able to produce explicit equations. In order to check his construction and improve the results, 
he made extensive use of SINGULAR to compute standard bases for global as well as for local 
orderings. One of his examples is the following: 

Example: ( Lossen, 1999 ) The irreducible curve with affine equation f(x,y) = 0, 



f(x,y) = y 



2y{x w • 



• r y 

21 



1 8 4 , 1 7 6 5 68, ' 5 10 

— x y H x y x y H x y 

8 y 16 y 128 y 256 y 



1024 



4 12 33 3 14 

x y + or ,^ Q x y 



429 2 16 

x y 

32768 y 



262144 

has degree 21 and an A228 _ singularity (x 2 



2048 

2U\ i 20 , iy 
y ) + x + x y 

— 0) as its only singularity. 



715 18 
xy 



65536 



228-smgLiiariiy [x — y 
In order to verify this, one may proceed, using SINGULAR, as follows 



ring s = 0,(x,y),ds; 

poly f = y2-2xl0y-x9y3+l/4x8y5-l/8x7y7+5/64x6y9-7/128x5yll+21/512x4yl3 

-33/1024x3yl5+429/16384x2yl7+x20-715/32768xyl9+xl9y2+2431/131072y21; 
matrix Hess = jacob(jacob(f )) ; //the Hessian matrix of f 

print (subst(subst (Hess ,x, 0) ,y,0) ) ; //the Hessian matrix for x=y=0 

==> 0,0, 
==> 0,2 

vdim(std( jacob(f ) ) ) ; //the Milnor number of f 

==> 228 



Since the rank of the Hessian at is 1, / has an Ak singularity at 0; it is an A22& singularity 
since the Milnor number is 228. To show that the projective curve C defined by / has no other 
singularities, we have to show that C has no further singularities in the affine part and no 
singularity at infinity. The second assertion is easy, the first follows from 

dimc(K[x,y]( X:V )/{}acob(f),f) = dim c (K[x,y]/{]acob(f), f), 

confirmed by SINGULAR: 

vdim(std(jacob(f )+f )) ; 

==> 228 //multiplicity of Sing(C) at (local ordering) 

ring r = 0,(x,y),dp; 
poly f = fetch(s,f); 
vdim(std(jacob(f )+f )) ; 

==> 228 //multiplicity of Sing(C) (global ordering) 



8. What else is needed 

In this survey I could only touch on a few topics where computer algebra has contributed to 
mathematical research. Many others have not been mentioned, although there exist powerful 
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algorithms and efficient i mplementations. 



for group action s of finite (Sturmfels, 1993 



(Derksen, 1997) or some uni-potent (Greue 



In the first pla c e, the computation of inva riant rings 
Kemper, 199£j; Decker and Dc Jong 1998 ), reductive 

prnnns hplonp' 



Ptister and "gchonemann, 1990-1998| " groups belong 
here. Computation of invariants have important applicati ons for explicit construction of modul i 
spaces, for example, for vector bundles or for singul arities ( Fruhbis-Kri iger, 2000; T. Bayer, 2000) 
but also for dynamical systems with symmetries (Gatermann, 1999). Libraries for computing 
invariants are available in SINGUL AR,. Available is also the Fuiseux expansion (even better, the 
Hamburger-Noether expansion, cf. Lamm (1999) for description of an implementation) of plane 
curve singularities. The latter is one of the few examples of an algorithm in algebraic geometry 
where Grobner bases are not needed. 

The applications of computer algebra and, in particular, of Grobner bases in projective alge- 
braic geometry are so numerous that I can only refer to the textbooks of Cox-Little-O'Shea, 
Eisenbud and Vasconcelos and the literature cited there. The applications include classifica- 
tion of varieties and vector bundles, cohomology, moduli spaces and fascinating problems in 
enumerative geometry. 

However, there are also some important problems for whi ch an algorithm i s either not known 
or not yet implemented (for further open problems see also Eisenbud (1993)[ ): 



Resolution of singularities 

This is one of the most important tools for treating singular varieties. At least three 
approaches seem to be possible: 

For surfaces we have Zariski's method of successive normalisation and blowing up points 
and the Hirzebruch-Jung method of resolving the discriminant curve of a projection. 
For arbitrary varieties, new methods of Bierstone, Milman and Villamajor provide a con- 
structive approach to resolution in the spirit of Hironaka. First attempts in this direction 
have been made by Schicho. 
Computation in power series rings 

This is a little vague since I do not mean to actually compute with infinite power series, 
the input should be polynomials. However, it would be highly desirable to make effective 
use of the Weierstrafi preparation theorem. This is related to the problem of elimination 
in power series rings. 

Moreover, no algorithm seems to be known to compute an algebraic representative of the 
semi-universal deformation of an isolated singularity (which is known to exist). 
Also, I do not know any algorithm for Hensel's lemma. 
Dependence of parameters 

In this category falls, at le ast principally, the study of Grobner bases over rings . 



of course, been studied, cf. 



Adams and Loustaunou (1994) , Kalkbrcnner (1998) 



This has, 
but I still 



consider the dependence of Grobner bases on parameters as an unsolved problem (in the 
sense of an intrinsic or predictable description, if it exists). 

In many cases, one is interested in finding equations for parameters describing precisely the 
locus where certain invariants jump. This is related to the above problem since Grobner 
bases usually only give a sufficient but not necessary answer. 

Mainly in practical applications of Grobner bases to "symbolic solving", parameters are 
real or complex numbers. It would then be important to know, for which range of the 
parameters the symbolic solution holds. 
Symbolic-numeric algorithms 

The big success of numerical computations in real life problems seems to show that sym- 
bolic computation is of little use for such problems. However, as is well-known, symbolic 
preprocessing of a system of polynomial (even ordinary and partial differential) equations 
may not only lead to much better conditions for the system to be solved numerically but 
even make numerical solving possible. 

There i s continuous progr ess in this direction, cf. Cox, Little and O'Shea (1998), JVIollei 



(1998), Verschelde (1999), not only by Grobner basis methods. A completely different ap 



proach via multivariate resultants (cf. Canny and Em iris (1997)) has become favourable 



to seve ral people due to the new sparse resultants by KJ clfand, Kapranov and Zelevinks 



(1994). However, an implementation in SINGULAR (cf"VVenk (1999; 7 Hllle brand (1999)) 



i 
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does not show superiority of resultant methods, at least for many variables against trian- 
gular set methods of either Lazard or Moller. Nevertheless, much more has to be done. 
The main disadvantage of symbolic methods in practical, real life applications is its com- 
plexity. Even if a system is able to return a symbolic answer in a short time, this answer 
is often not humanly interpretable. Therefore, a symbolic simplification is necessary, ei- 
ther before, during, or after generation. Of course, the result must still be approximately 
correct. 

This leads to the problem of validity of "simplified" symbolic computation. A completely 
open subproblem is the validity resp. error estimation of Grobner basis computations with 
floating point coefficients. 

The simplification problem means providing simple and humanly understandable symbolic 
solutions which are approximately correct for numerical values in a region which can be 
specified. This problem belongs, in my opinion, perhaps to the most important ones in 
connection with applications of computer algebra to industrial and economical problems. 
Non-commutative algorithms 

Be fore Grobner bases were introduced by Buchberger, the so-called Ritt-Wu method, 
cf. |Ritt (1950) , was developed for symbolic computation in non-commutative rings of 
differential operators. However, nowadays, commutative Grobner bases are implemented 
in almost every major computer algebra system, whilst only few systems provide non- 
commutative algorithms. Standard bases for some non-com mutative structures have been 
implemented in the system FELIX (Apel and Klaus, 1991) as well as in an experimental 
version in SINGULAR; the system Bergman and an extension called Anick can compute 
Grobner bases and higher syzygies in the non-commutative case. 

Highly desirable are effective implementations for non-commutative Grobner bases in the 
Weyl algebra, the Grassmannian, for D-modules or the enveloping algeb ra of a finite 



dimensional Lie al gebra ("the ge neral theory being basically understood, cf. Mora (1989) 
Ufnarovski (1998)[ jApel (1998)|). 

The recent textbook of [jaito, Sturmfels and Takayama (1999) shows a wide variety of 
algorithms for modules over the Weyl algebra and W-modules for which an efficient im- 
plementation is missing. 

But ev en classical algebraic geometry, as was shown, for example, by Kontsevich and Manin 



(1994), has a natural embedding into non-commutative algebraic geometry. A special 



case is known as quantisation, a kind of non-commutative deformation of a commutative 
algebra. 

Providing algorithms and implementations for the use of computer algebra in non-com- 
mutative algebraic geometry could become a task and challenge for a new generation of 
computer algebra systems. 
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